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The confinement improvement in reversed-shear experiments on the Tokamak Fusion Test Reactor 
[Plasma Phys. Controlled Fusion 26, 11 (1984)] is investigated using nonlinear gyrofluid 
simulations including a bounce-averaged trapped electron fluid model. This model includes 
important kinetic effects for both ions and electrons, and agrees well with linear kinetic theory. Both 
reversed shear and the Shafranov shift reverse the precession drifts of a large fraction of the trapped 
electrons, which significantly reduces the growth rate of the trapped electron mode, found to be the 
dominant instability in the core. Two positive feedback transition mechanisms for the sudden 
improvement in core confinement are discussed: (1) Shafranov shift suppression of the trapped 
electron mode, and (2) turbulence suppression by radially sheared ExB flows. While both effects 
appear to be playing roles in the transition dynamics in most experiments, we show that Shafranov 
shift stabilization alone can cause a transition. © 1997 American Institute of Physics. 
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I. INTRODUCTION 

In this paper, we investigate the improved confinement 
with reversed magnetic shear 1 ' 2 using nonlinear gyrofluid 
simulations extended to higher accuracy 3 and including a 
bounce-averaged trapped electron fluid model. 4 This nonlin¬ 
ear electron fluid model includes the kinetic effects of the 
trapped-electron precession resonance and retains the full 
pitch-angle dependence of the electron response. Retaining 
the pitch-angle dependence is essential to describe the sup¬ 
pression of trapped electron modes (TEM) in enhanced 
reversed-shear (ERS) discharges, where the dominant stabi¬ 
lizing effect is the reversal of the toroidal precession drifts of 
barely trapped electrons. This model has been validated by 
detailed linear comparisons with the most comprehensive ki¬ 
netic calculations, 5-7 and has extended our simulations to 
include the ERS and supershot core in the Tokamak Fusion 
Test Reactor (TFTR), 8 where the dominant instability is the 
TEM. In the core of reversed-shear RS plasmas (before the 
transition to ERS) and supershots, our simulations predict 
fluxes in rough agreement with TRANSP, 9 and the ion heat 
transport is convection dominated. After the transition to 
ERS, the TEM is strongly suppressed by the combination of 
negative magnetic shear and the Shafranov shift. In ERS, the 
Shafranov shift becomes very large, and is more important 
than negative magnetic shear in reversing the drifts of the 
trapped electrons. After the transition, a shorter wavelength 
TEM is still weakly unstable, but nonlinear simulations find 
that the transport is reduced by about a factor of 40, in rough 
agreement with experiment, although the electron heat trans¬ 
port is underestimated. In contrast to magnetic shear, the 
Shafranov shift stabilization is a positive feedback mecha¬ 
nism and is a possible trigger for the sudden transition to 
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ERS, since steeper pressure gradients lead to larger shifts, 
more drift reversal, less transport, and in turn, to even steeper 
pressure gradients. Another potentially important mechanism 
is stabilization via radially sheared electric fields. 10-13 We 
find that in ERS, the amount of electric-field shear is usually 
near the level required to completely stabilize the TEM, but 
that without the Shafranov shift stabilization, the growth 
rates are too high for the measured levels of ExB shear to 
stabilize the turbulence. Therefore, both stabilization mecha¬ 
nisms seem to be important in causing the transition. Further 
studies are needed to elucidate the relative importance of 
these two mechanisms. 

II. BASIC MODEL AND LINEAR RESULTS 

To accurately describe the core confinement improve¬ 
ment in ERS discharges, kinetic effects are crucial. Our basic 
model consists of six-moment toroidal ion gyrofluid 
equations 3 for each ion species (typically D, C, and a Max¬ 
wellian beam component), and bounce-averaged fluid equa¬ 
tions for the trapped electrons, 4 coupled by quasineutrality. 
These equations employ the electrostatic approximation, 
which is supported by previous results 6 ' 7 showing that the 
dominant finite-/? effect comes through variations in the 
magnetohydrodynamic (MHD) equilibrium, and the effects 
of coupling to magnetic fluctuations are much smaller for the 
/? values considered here. The toroidal ion equations include 
the kinetic effects of parallel 14 and toroidal drift phase mix¬ 
ing, finite Larmor radius effects, 15 and trapped ion effects. 
The electron equations include toroidal precession drift 
phase mixing, and because the electron moments are func¬ 
tions of pitch angle, we use a pitch-angle scattering operator 
for electron-electron and electron-ion collisions. 4 With this 
collision operator, the electron equations are continuously 
valid from the collisionless regime, where the nonadiabatic 
electron response is driven by the trapped electron precession 
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FIG. 1. Comparison of linear growth rates from the gyrofluid model (solid) 
and fully kinetic results of Ref. 5 (dashed) as collisionality is varied, dem¬ 
onstrating good agreement. 


over k 8 pi. In the core region, r/a <0.45, the dominant insta¬ 
bility is a high-A' e trapped electron mode, and a lower k g ion 
temperature gradient (ITG) driven mode is dominant for 
rla >0.45. When the trapped electrons are turned off by as¬ 
suming an adiabatic electron response, the measured profiles 
are stable for rla <0.45. This behavior is typical of RS 
modes, ERS modes, and supershots: the TEM is the domi¬ 
nant instability in the core region where the density gradients 
are steep, outside this region the ITG mode dominates, and 
also in low-confinement modes (L modes). Both the kinetic 
and gyrofluid calculations use general magnetic geometry, 
but the kinetic calculation uses a slowing-down beam distri¬ 
bution function while the gyrofluid calculation uses a Max¬ 
wellian beam distribution. The ITG mode is stable in the 
core because of the steep V n , and low 77 , . Although the 
TEM in the core has a larger peak growth rate, it causes 
much less transport than the ITG mode, which is unstable for 
rla >0.45, since the ITG mode is unstable at significantly 
longer wavelengths. Across the transition from RS to ERS, 
the longest wavelength trapped electron modes are stabi¬ 
lized, reducing the transport. 


resonance, to the plateau collisional regime where the elec¬ 
trons become adiabatic. These equations have been carefully 
benchmarked against linear kinetic theory. Figure 1 shows a 
comparison of the linear growth rate from this model with 
fully kinetic calculations 5 for the parameters 77 ^= 77 , =3, 
L„/R = 1/3, e = 1/6, s = l, <7 = 1.5, and /3 = 0, as the collision¬ 
ality is varied, demonstrating good agreement. 

We now examine the linear instabilities in an ERS mode. 
Figure 2 shows linear growth rates and real frequencies ver¬ 
sus minor radius well after the transition to ERS. The mea¬ 
sured experimental parameters at each radius (from TRANSP) 
are used as inputs for the calculations. The eigenfrequencies 
shown are for the fastest growing mode, with y maximized 
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FIG. 2. Maximum linear growth rates and corresponding real frequencies vs 
r /a for an ERS discharge well after the transition, again demonstrating good 
agreement between the gyrofluid model (solid) and fully kinetic calculations 
(dashed). 


III. FINITE-/? STABILIZATION OF THE TRAPPED 
ELECTRON MODE 

We now discuss the physical mechanisms that stabilize 
the trapped electron mode in the core. Both negative mag¬ 
netic shear and the Shafranov shift cause favorable preces¬ 
sion of all but deeply trapped electrons and can stabilize the 
TEM. The precession reversal due to negative s is relatively 
well known, and its stabilizing influence has been thoroughly 
discussed recently . 16 The finite-/? drift reversal effect was 
first pointed out in Ref. 17. We begin by using the s — a 
model equilibrium to demonstrate this effect, where the 
Shafranov shift effects are measured by a =—q 2 Rd/3ldr. 
The precession drift is the bounce average of the VZ? and 
curvature drifts, and in the s — a model is given by 

k f)P pX) fp 

( w de)o( K ) = — n —( cos 0+(sff—a sin #)sin 9) b 

K 


= eP ft ' e [G 0 (k) + sG s (k) +aG a (K)\, (1) 

where G q (k)=2E(k 2 )SK(k 2 )-1, G s (k)=4[E( k 2 )IK( k 2 ) 
+ k 2 - 1 ], andG a (/c) = (4/3)[(l - 2k 2 )E(k 2 )/K(k 2 ) + k 2 
— 1 ], where E and K are complete elliptic integrals, and k is 
a pitch-angle variable which is zero for deeply trapped elec¬ 
trons and one for barely trapped electrons, as defined in Ref. 
4. Figure 3 shows these three terms. Since G s > 0 and G a < 
0, both negative s and positive a reverse the precession of 
barely trapped electrons. This drift reversal is stabilizing, as 
the reversed electrons can no longer resonate with the TEM. 

For TFTR ERS core parameters, the Shafranov shift in¬ 
duced drift reversal (a) actually dominates, since the shear is 

only weakly negative, s -1/4, but the shift becomes very 

large, a~ 2, typical parameters of the second stability re¬ 
gime to MHD ballooning modes. We show this in Fig. 4 by 
varying s and a , keeping other parameters fixed, for 
rla =0.25 at t = 3 s of TFTR shot 84011. Using s = 0.25 and 
a = 0.3, values typical of the supershot regime, very few 
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FIG. 3. Pitch-angle dependence of the three terms in the s — a precession 
frequency. Both negative shear (s<0) and the Shafranov shift (a>0) cause 
favorable precession of barely trapped electrons. This is stabilizing, as these 
drift reversed electrons can no longer resonate with the trapped electron 
mode. 


electrons have reversed precession, as shown in Fig. 4(a). 
The corresponding growth rates for the trapped electron 
mode are shown in Fig. 4(b). Reversing the magnetic shear 
to values typical in the reversed-shear mode before the tran¬ 
sition (RS), s = —0.25, but before the shift has become large, 
a = 0.3, reverses the precession drifts of more of the barely 
trapped electrons and stabilizes the TEM somewhat. Using 
the measured parameters after the transition, s= — 0.25 and 
a = 2.5, we find much more drift reversal. The long wave¬ 
lengths are stabilized, and the remaining high-A TEM is quite 
weak and causes little transport. Thus, the effect on the 
growth rates of reversed shear alone is rather weak, but the 
combined effects of reversed shear and the Shafranov shift 
are very strongly stabilizing. 

To more accurately describe this effect, we have ex¬ 
tended our simulations to general magnetic geometry. Using 
the output of an equilibrium code, we numerically compute 
the dependence along the field line coordinate, ft of the geo¬ 
metric terms that enter the gyrofluid equations, as described 
in Ref. 18, u> d (9), A ± (ft), b-V( ft), and B(9). In addition, 
we numerically bounce average co d ( 9) to calculate the toroi¬ 
dal precession frequency for the trapped electrons: 
(w d ) b (K). To show the effects of full geometry on the pre¬ 
cession frequency, we calculate equilibria using JSOLVER, 19 
and the measured TFTR profiles after the transition, at 
t = 3.0 s. We can also artificially reduce the Shafranov shift, 
by repeating the equilibrium calculation with all densities 
reduced by a factor of 10. We then numerically calculate the 
bounce-averaged curvature and VB drifts using the JSOLVER 
output. The toroidal precession frequencies are compared for 
these two cases in Fig. 5. Because the q profile changes 
slightly when the density is reduced, we compare at rla 
= 0.3, where s= — 0.17 and A' = 0.53 for the full density 


(a) K 



FIG. 4. Precession frequency (a) and corresponding growth rates vs k Q (b) 
as s and a are varied. For typical supershot values (s = 0.25 and a = 0.3), 
very few trapped electrons precess in the favorable direction and the TEM 
growth rate is fairly strong. For RS values, where the shear is reversed but 
the shift is not very large (s= —0.25 and a = 0.3), more electrons precess in 
the favorable direction, and the TEM growth rate is reduced somewhat. For 
ERS values, where the shift has become quite large (s=— 0.25 and 
ai = 2.5), the majority of trapped electrons precess in the favorable direction, 
and the TEM growth rates are dramatically reduced. The longer wave¬ 
lengths, which cause more transport, are completely stabilized. 

case (A' = dA/dr, where A is the Shafranov shift), and at 
rla = 0.26 for the reduced density case, where the shear is 
the same, s= —0.17, and A' =0.07. The effect of the Shafra¬ 
nov shift is quite striking, and is significantly different from 
that given by the s — a model, though the basic trend is simi¬ 
lar. Because of the pitch-angle (/c) dependence of the Jaco¬ 
bian of the transformation to v , k variables, the fraction of 
electrons with drift reversal is given by n e /n e 0 
= 1 - f - e (d9/4) Vl - (B/B mm )(\ - 2e B Ko), where ( w de)b( K o) 
= 0, 9, is the turning point for a particle with k= k 0 , and the 
other notation is defined in Ref. 4. From Fig. 5, we find 75% 
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FIG. 5. Toroidal precession frequency vs pitch angle using full geometry, 
for A ' = 0.53 (solid) and 0.07 (dashed). In full geometry, the Shafranov shift 
is dramatically stabilizing. In contrast to the s — a model, the precession 
frequency of all trapped electrons is reduced by A', for all pitch angles. 


drift reversal at A' = 0.53 and 60% drift reversal at 
A' =0.07. We have also calculated the case with A' = 0.07 
and s = 0, at r/a = 0.29, and find 40% drift reversal. 

Some of the differences between the full geometry re¬ 
sults and the s — a model can understood by looking at the 
precession drift frequency in the low-/3, small e = r/R limit, 
where (3~0(e 2 ) and A '~0(e). The combined curvature 
and V B drift frequency is (for 9 0 = 0) 

2 

iw d $> = 2 ^-gZ [BxV S + 6 2 bx ( b Vb )] -V(2) 

and for low (3 and small e, this becomes 

e 

cos 9— (A' + e)H—i+sZZsin rA" siir 9 

+ A's9 sin 0 cos 6— s(A' + e)sin 2 9+0(e 2 ). 

(3) 

Comparison with Eq. (1) shows that the s— a model only 
keeps three of these terms, the cos 0 term, the s 9 sin 8 term, 
and the — rA" sin 2 0 term, which corresponds to the 
a sin 2 9 term in Eq. (1), because A' °c (Z,/2 + /3 g )r and 
a= — q 2 R(d /3/dr) ' y A". The s — a model misses the reduc¬ 
tion in precession frequency for all particles (not just the 
deeply trapped particles) from A'. This A' reduction of a> d 
arises from compression of the flux surfaces and a change in 
the tilt of the field lines, which comes from variations in 
B p over the surface. This reduces k g in the bad curvature 
region and increases k g in the good curvature region, in ad¬ 
dition to a sin 0 variation in k r proportional to A'. These 
pieces combine with the 9 cos 9+v sin 9 variation of 
BxVB to give the constant A' reduction in u> d (9), which 
after bounce averaging leads to a constant A' reduction in 


( a> d ) b (K ), independent of pitch angle. In this ordering, 
f3~rt{e 2 ), which is reasonably well satisfied for TFTR, the 
changes in B are subdominant. The Shafranov shift also re¬ 
duces the field line length in the bad curvature region and 
increases the field length in the good curvature region, but 
this does not affect the constant A' reduction in (oj d ) h , since 
the bounce-averaged drift {<a d ) b = (^JZw d /|i;|||)/(^JZ/|i;|||). 
Using dl = {rB/Bp)d8 qRo(l — A'cos 0)d9, we find r b 
= §dll\v^\\-(A'-€/2 )Gq\K(k 2 ), and 

(^) fc =Go-(A' + e )+A + (A'-|)(G 2 -l-G a ) 

+ sG s j 1 + A'G 0 —| G 0 J +(rA" + sA' +se)G a 

+ ^(2G S (1—2/r 2 ) —GJ + G(e 2 ), (4) 

where the Gs are defined following Eq. (1) and are plotted in 
Fig. 3. The unusual feature of the Shafranov shift stabiliza¬ 
tion of the TEM is that it becomes more stabilizing for 
steeper pressure gradients, and allows access to what could 
appropriately be called the second stability regime to the 
TEM. This feature also makes Shafranov shift stabilization a 
potential positive feedback mechanism in producing the sud¬ 
den transition to improved confinement in ERS. 

IV. ERS TRANSITION MECHANISMS 

To get a sudden transition or bifurcation to the improved 
confinement in ERS, transport must decrease as the driving 
gradients are increased, providing the positive feedback nec¬ 
essary for the core density to run away. There appear to be 
two potential mechanisms that can produce this positive 
feedback: (1) The stabilization of the TEM from the Shafra¬ 
nov shift, 20 and (2) radial ExB flow shear stabilization of 
the turbulence. 10 As discussed below, it appears that both of 
these effects may produce a transition individually, or in tan¬ 
dem. Both of these effects are stronger in ERS plasmas than 
in other tokamak regimes. ERS modes have unusually large 
A' since both q and (3 are large in core, so Shafranov shift 
stabilization is enhanced. This becomes a positive feedback 
mechanism when the Shafranov shift stabilization is strong 
enough to overcome the increased instability drive from the 
steeper gradients. Then as V p increases. A' increases, the 
transport decreases, and the pressure gradients steepen fur¬ 
ther. 

When the shearing rate, a > E , from radial ExB flow 
shear is on the order of or larger than the maximum linear 
growth rate, w E >y, the turbulence can be completely 
stabilized. 12 When the radial electric field is driven by the 
pressure gradient, this is also a positive feedback mecha¬ 
nism, since as V p increases, E' r increases, the transport de¬ 
creases, and the pressure gradients steepen further. This ef¬ 
fect is also strong in the ERS regime, since ERS modes have 
large Vp. The appropriate shearing rate in general 
geometry 11 is u> E =(RB p /B)d/dr(E r /B p R), and since B p is 
small in the core of ERS modes, the shearing rate is en¬ 
hanced. In addition, the Shafranov shift enhances u> E . Tur- 


Phys. Plasmas, Vol. 4, No. 5, May 1997 


Beer et al. 


1795 






FIG. 6. Comparison of maximum linear growth rates (closed circles) and 
measured ExB shearing rate, <o E (open circles). Before the transition, the 
maximum linear growth rate exceeds the ExB shearing rate. After the tran¬ 
sition, the linear growth rate is significantly less than the E X B shearing rate, 
so the turbulence is suppressed. Without the stabilizing effect of A' 
(dashed), however, the linear growth rate increases as the density gradient 
steepens, staying above a> E . Thus, in this case, both effects appear to be 
playing roles in the transition dynamics. 

bulence suppression by radial ExB flow shear has been 
shown to be consistent with the transition to improved con¬ 
finement in ERS in Ref. 21. 

Figure 6 compares the measured ExB shearing rate 
within the good confinement zone, at r = 0.25a, to the maxi¬ 
mum linear growth rate calculated at the same radius, for an 
ERS mode that transitioned at t = 2.71 s. The measured 
ExB shearing rate is larger than y after the transition, dem¬ 
onstrating that a> E is large enough to substantially suppress 
the turbulence. Also shown are the linear growth rates with¬ 
out Shafranov shift stabilization (A'=0), where the linear 
growth rates are too large for a> E to completely stabilize the 
turbulence. Therefore, in this case, it appears that both 
mechanisms are playing a role: Shafranov shift stabilization 
keeps the linear growth rates from increasing with the pres¬ 
sure gradient, and eventually u> E overcomes y lin . It is diffi¬ 
cult to perform experiments that clearly separate these two 
effects, since both are proportional to V p, and both A' and 
E\. usually increase dramatically across the transition. Back- 
transition experiments on TFTR, 22 varying the co-counter 
beam fraction, show that decreasing E' r can cause the plasma 
to drop out of ERS, indicating that E' r is playing a key role in 
maintaining the transport barrier. There are cases, however, 
where the criterion u> E > y does not correlate well with the 
transition. 23 The observed B scaling of the threshold power, 
ocR 2 , also seems more consistent with a transition induced 
by Shafranov shift stabilization, since A'°c/3. 

We now describe a theoretical experiment, which shows 
that the Shafranov shift appears to be able to produce a tran¬ 
sition independently, without a> E stabilization. To cause a 
transition, Shafranov shift stabilization must overcome the 




FIG. 7. Calculated particle flux as Vzz is varied, (a) Within the good con¬ 
finement zone, at r/a = 0.3. For low-density gradients transport increases 
with the gradient, but beyond a critical density gradient of about 0.6 times 
the measured value at f=3 s, Shafranov shift induced stabilization of the 
trapped electron mode overcomes the increasing density gradient, and the 
particle flux decreases. Thus, beyond 0.6, the Shafranov shift becomes posi¬ 
tive feedback, causing a transition. Once this critical gradient is exceeded, 
the density would run away, (b) Outside the good confinement zone, at 
rla = 0.5, the particle flux is still increasing with increasing density gradi¬ 
ents so there is no transition. This correlates well with experiment, where 
the density is increasing linearly in time only within r/a<0.4. 


increased transport from steeper gradients. To test this effect, 
we generate a series of equilibria using JSOLVER, 19 starting 
with the measured profiles at t = 3 s, and gradually decreas¬ 
ing all densities. This gives a series of equilibria with A' 
r s- n, and very similar q profiles. We then calculate the par¬ 
ticle flux, T= —DVn, for each of these equilibria, either by 
using fully nonlinear simulations or by estimating the par¬ 
ticle diffusivity D from the maximum over k 8 of y/k 2 , 
which correlates well with our nonlinear simulations. 24 The 
results are shown in Fig. 7, using D = (5/3) y/k from linear 
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FIG. 8. Calculated particle flux as Vn is varied, now only changing the 
densities within the good confinement zone, r/a< 0.4. Again, beyond some 
critical density gradient, Shafranov shift stabilization of the TEM overcomes 
the increasing instability drive from the increasing gradient, causing a tran¬ 
sition. Including co E stabilization would decrease the critical gradient, as 
shown schematically (dashed line), by shifting the peak of the flux vs V/? 
curve. 


calculations with the s —a model and with full geometry, 
and two nonlinear simulations in full geometry. The two 
nonlinear simulations agree well with the y/k ± estimate with 
a proportionality constant of 5/3. Figure 7(a) is within the 
good confinement zone, at r/a = 0.3. At small density gradi¬ 
ents, we find the usual situation where steeper density gradi¬ 
ents cause more transport. But at roughly 60% of the mea¬ 
sured density gradient at t = 3 s, the stabilization from the 
Shafranov shift overcomes the increasing density gradient 
and causes a transition, since beyond this point steeper den¬ 
sities cause less transport and the density runs away. Figure 7 
shows the same calculation, at rla = 0.5, outside the good 
confinement zone. Here, transport continues to increase with 
increasing density, so no runaway occurs. Thus, this result 
correlates well with the observed transport barrier location in 
the experiment. 

In the experiments, however, the entire density profiles 
do not increase; the density only increases in the core. We, 
thus, repeat the above experiment only changing the density 
in the good confinement zone, out to the “knee” in the den¬ 
sity profile at r/a = 0.45. This series also decreases 
L n = \(l/n 0 )(dn 0 /dr)]- 1 , and increases the instability drive 
at the same time that A' is increased. These results are shown 
in Fig. 8 using full geometry. The Shafranov shift still over¬ 
comes the increased drive and causes a transition at about 
half the measured gradient. In Fig. 8 we have also schemati¬ 
cally indicated the effects of u > E , which are not included in 
these calculations, by reducing the flux by a contribution 
proportional to V«/(l— A') 2 , including the Shafranov shift 
enhancement of Since ExB shear is also °c V n, the 
stabilization from u> E increases with V n, and reduces the 



FIG. 9. Nonlinear simulations of an ERS discharge at rla = 0.25 (a) before 
the transition, and (b) after the transition. Before the transition, the predicted 
fluxes are in reasonable agreement with power balance measurements, and 
the ion heat flux is convection dominated. After the transition, the transport 
is greatly reduced, as the longest wavelengths are stabilized by the large 
Shafranov shift. 


threshold gradient by shifting the peak of the flux versus 
V« curve to lower V«. 

V. NONLINEAR SIMULATIONS 

We now describe nonlinear toroidal gyrofluid simulation 
results, including deuterium, carbon, and trapped electrons 
for an ERS discharge within the good confinement zone, at 
rla = 0.25, where the turbulence is driven by trapped elec¬ 
tron modes. The results just before the transition, at 
t = 2.66 s, are shown in Fig. 9(a). We find reasonable agree¬ 
ment with measured fluxes: the predicted fluxes are about 
three times the measured values from TRANSP. In addition, 
the ion heat flux is convection dominated, as seen experi- 
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mentally. A simulation using measured input parameters af¬ 
ter the transition is shown in Fig. 9(b). The turbulent fluxes 
drop by about a factor of 40 compared to the pretransition 
simulation, and the fluctuation levels drop dramatically. 
However, this run has not yet reached a satisfactory steady 
state, which would require a much longer run with these 
small growth rates. Another caveat is that so close to mar¬ 
ginal stability, the simulation results are very sensitive to 
input gradients. Although the electron heat flux is below ex¬ 
perimental measurements after the transition, the dramatic 
confinement improvement across the transition is repro¬ 
duced. 

VI. CONCLUSIONS AND DISCUSSION 

To summarize, we have investigated the confinement 
improvement in ERS regimes using nonlinear gyrofluid 
simulations. Our comprehensive gyrofluid model 3 ' 4 agrees 
with kinetic calculations, where we find that a high-/r w 
trapped electron mode is the dominant instability in the steep 
density gradient core region, and that a lower-/: # ITG mode 
is unstable for r/a> 0.45. Our pitch-angle dependent trapped 
electron model was required to capture precession drift re¬ 
versal, which greatly reduces the growth rate of the TEM in 
the core. The Shafranov shift accounts for much of this drift 
reversal, and we have extended our simulations to general 
geometry to accurately capture this effect. Shafranov shift 
stabilization is also a potential positive feedback mechanism 
to produce the sudden transition to enhanced confinement in 
ERS. We also show that turbulence suppression via radial 
electric-field shear is often important, and usually find 
Tim~ (j) k at the transition. It appears that both effects are 
working together to produce the transition: Shafranov shift 
stabilization keeps the growth rates from increasing with the 
driving gradients, and keeps the turbulence weak enough to 
allow the increasing pressure driven ExB shear to further 
stabilize the turbulence. Finally, our nonlinear simulations 
capture the main features of the transition to ERS, even 
though they do not include equilibrium scale sheared ExB 
flows. The calculated pretransition fluxes are in rough agree¬ 
ment with TRANSP, the core ion heat transport mostly con¬ 
vective, and we find a sharp drop in transport at transition, as 
the long wavelengths are stabilized, and the remaining high- 
k TEM causes little transport. 

Many questions for future work remain. Perhaps some 
high-/: mode or subcritical dco pe scale turbulence is causing 
the residual electron transport after the large scale turbulence 
is stabilized. We plan to include E' r in our nonlinear gyro¬ 
fluid simulations to investigate parametric dependences of 
the approximate criterion y lin ~ a > E . A crucial question for 
future experiments is whether we can find a regime where 
Shafranov shift stabilization alone gives turbulence suppres¬ 
sion. The Shafranov shift has a more favorable scaling to 
reactors than E \., since the Shafranov shift is independent of 
machine size, while E\. effects become weaker in larger ma¬ 
chines. An exciting alternative is to drive E\. with ion Bern¬ 
stein waves instead of pressure or velocity gradients. We are 
also beginning to investigate other geometries that may tran¬ 
sition more easily. Oblate reversed triangularity 25 tokamaks. 


which have stronger Shafranov shifts at the same (3 as con¬ 
ventional tokamaks, are one possibility. Low aspect ratio to¬ 
kamaks, such as the proposed National Spherical Torus Ex¬ 
periment, have naturally large drift reversal from low R/a 
and high A', and Shafranov shift stabilization should be quite 
strong. Finally, perhaps stellarator configurations can be op¬ 
timized for drift reversal to take advantage of the confine¬ 
ment improvement in ERS. Our results indicate that obtain¬ 
ing steep density gradients may be crucial, since it appears to 
be easier to stabilize the TEM (through the Shafranov shift 
and precession drift reversal) than the ITG mode. 
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